Supergalactic winds driven by multiple superstar clusters. 
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ABSTRACT 

Here we present two dimensional hydrodynamic calculations of free expanding supergalactic 
winds, taking into consideration strong radiative cooling. Our main premise is that supergalac- 
tic winds are powered by collections of superstar clusters, each of which is a source of a high 
metallicity supersonic diverging outflow. The interaction of winds from neighboring superstar 
clusters is here shown to lead to a collection of stationary oblique shocks and crossing shocks, 
able to structure the general outflow into a network of dense and cold, kpc long filaments that 
originate near the base of the outflow. The shocks also lead to extended regions of diffuse soft 
X-ray emission and furthermore, to channel the outflow with a high degree of coUimation into 
the intergalactic medium. 
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1. Introduction 

Supergalactic winds (SGWs; see Heckman 
2001), as first envisaged by Chevaher and Clegg 
(1985; hereafter CC85), freely flow from the nu- 
clear starburst region of a galaxy into the inter- 
galactic space. The flow originally thought to 
behave in an adiabatic manner has recently been 
shown to be strongly affected by radiative cool- 
ing, particularly for powerful (> lO''^ erg s~^) and 
compact starbursts (see Sihch et al. 2003; here- 
after referred to as Paper I). Radiative cooling 
does not disturb the velocity of the outflow, nor 
the expected density drop (~ r^^) acquired as the 
flow moves away from the massive central cluster. 
Rapid cooling however, brings the temperature 
down in regions close to the starburst, favoring 
rapid recombination, and thus making the fast 
streaming outflow easy target of the UV radiation 
generated by the central stars. The prediction is 
then a much reduced size of the X-ray emitting 
zone and a fast moving HII region gas originating 
close to the central starburst. 

M82 with its wind striking the "Ha cap" at 
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10 kpc from the central source (see Devine & 
Bally 1999), is without a doubt the best exam- 
ple of a SGW in the Local Universe, injecting 
its newly processed metals into the intergalactic 
space (IGM). The central biconical wind in M82 
however, presents an intricate structure that has 
little to do with the outcome from models of su- 
perwinds. In particular the X-ray and the Ha fil- 
amentary structure as well as their spatial coin- 
cidence have presently no explanation. As stated 
by Strickland & Stevens (2000) the predicted X- 
ray luminosity, even using the adiabatic model 
of CC85, falls short by more than three orders 
of magnitude below the observed value. Strick- 
land et al. (1997) also showed that the entropy 
of the X-ray emitting gas increases with distance 
from the plane of the galaxy, fact that is incon- 
sistent with an adiabatic outflow model. On the 
other hand, the Ha filamentary structure, beau- 
tifully evidenced by Subaru (see Ohyama et al. 
2002 and references therein), is clearly not limb 
brightened and originates right at the base of the 
outflow reaching several kpc in a direction almost 
perpendicular to the galaxy plane. These facts 
are also very different from the superbubble fea- 
tures calculated by Suchkov et al. (1994) in which 
a complex filamentary structure develops at large 
distances (several kpc) from the galaxy disc, as a 
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result of matter entraining the hot superbubble. 
There should be in fact very little resemblance be- 
tween supergalactic winds freely expanding into 
the IGM and the calculations mentioned above 
in which a supershell contains always the possi- 
ble outflow of the superbubble interior, and thus 
inhibits the development of a super galactic wind. 

So far, all calculations in the literature have as- 
sumed that the energy deposition arises from a 
single central cluster that spans several tens of pc, 
the typical size of a starburst. However, recent 
optical, radio continuum, IR and UV observations 
(Ho 1997; Johnson et al. 2001; Colina ct al. 2002; 
Larsen & Richtler 2000, Kobulnicky & Johnson 
1999) have revealed a number of unusually com- 
pact young stellar clusters. These overwhelmingly 
luminous concentrations of stars present a typical 
half-light radius of about 3 pc, and a mass that 
ranges from several times 10'' M0 to a few 10^ 
Mq. Clearly these units of star formation (super- 
star clusters, SSCs) are very different to what was 
usually assumed to be a typical starburst. Several 
of these entities have now been identified within 
a single starburst nucleus. For example, there are 
about 100 of them composing a flattened distribu- 
tion of 150 pc radius at the center of M82, (see de 
Grijs et al. 2001, O'Connell et al. 1995), at least 
four within the nuclear zone of NGC 253 (Watson 
et al. 1996), and many of them in the Antenae 
(Whitmore et al. 1999). 

Here we investigate the effects that the pres- 
ence of several of these young compact clusters in a 
galaxy nucleus may have on the inner structure of 
well developed, or freely expanding, SGWs. Sev- 
eral aspects are considered in this two dimensional, 
first approach, to the 3-D interaction of multiple 
powerful winds. Among these, the metallicity of 
the superwind matter is here shown to have a pro- 
found impact on the inner structure of SGWs. 

Here we study the interaction of several strong 
winds emanating from a collection of nearby su- 
perstar clusters, causing together the development 
of an extended region in which a plethora of cross- 
ing shocks collimate the general outflow while giv- 
ing rise to an important soft X-ray contribution 
at large distances from the starburst nuclei. An 
extended region in which the layers of strong di- 
rect wind interactions lead, under strong radiative 
cooling, to a well developed network of elongated 
filaments. Section 2 displays our two dimensional 



calculations that use CC85 as initial condition in 
each of the superstar clusters. The evolution leads 
to multiple interactions, and the outcome of these, 
depending on the local values of density, temper- 
ature and metallicity, define when the flow is af- 
fected by strong radiative cooling. Section 3 dis- 
cusses some of the observational consequences of 
such well coUimated and structured outflows, in 
particular the resultant filling factor for different 
gas phases. 

2. The interaction of winds from superstar 
clusters 

There are two different evolutionary stages in 
the remnants produced by the large mechanical 
energy input rate associated to nuclear starbursts. 
The first one, the superbubble phase, is that dur- 
ing which a large-scale remnant evolves within the 
ISM, either within the disk or within the halo of 
the host galaxy. The second phase, the freely 
expanding supergalactic wind, in which an open 
channel in the ISM allows the supernova matter 
to freely stream into the intergalactic medium. 

Weather one assumes a central source or multi- 
ple SSCs composing a starburst nucleus, the initial 
interaction with the host galaxy ISM will be rather 
similar. In particiflar, both sources of energy will 
lead to the development of a leading shock, able to 
sweep the surrounding ISM into a large-scale su- 
pershell, while a reverse shock will cause the full 
thermalization of the matter violently ejected by 
the central source or the multiple stellar clusters 
(see for example figure 1 in Heckman et al., 1990). 

In both cases the thermalized hot gas within 
the superbubble will power the leading shock and 
thus promote the growth of a remnant able to 
eventually exceed the dimensions of the galaxy 
disk, causing via Rayleigh - Taylor instabilities the 
fragmentation and disruption of the supershell of 
swept up matter (see figure 4 in Tenorio-Tagle & 
Bodenheimer, 1988). 

If one considers the presence of an extended 
halo (see Silich & tenorio-Tagle, 1998, 2001), the 
thermalized wind energy will generate an even 
larger supershell (see for example figure G in 
Strickland & Stevens, 2000), the velocity of which 
will continuously decay as more halo matter is in- 
corporated. The presence of extended haloes (see 
e.g. Melo et al., 2002) is also supported by the ex- 
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istence of large-scale supershclls in a large variety 
of galaxies (see for example Oey et al. 2002 and 
Marlowe et al., 1995). 

Note that a supergalactic wind does not develop 
until the remnant exceeds the dimensions of the 
galaxy disk, or in case of an extended halo, until 
the superbubble reaches the outskirts of the galaxy 
(for a review see Heckman, 2001). At this moment, 
the leading and the reverse shocks cease to exist 
and the metals ejected by the star clusters will 
freely move into the intergalactic medium causing 
its contamination. It is also at this stage, once a 
channel has been carved in the ISM, once the su- 
perwind has developed, that major differences will 
arise from the assumption of a single or a multiple 
source of energy. 

• In the case of a single compact source of en- 
ergy, the free expanding wind has recently 
been found to be subjected to strong radia- 
tive cooling (see Paper I). Thus, the adi- 
abatic model of Chcvahor & Clegg (1985), 
used to calculate the extended X-ray emis- 
sion of superwinds, leads to large overesti- 
mates. More realistic calculations, account- 
ing for radiative cooling, have shown a much 
reduced volume, of only a few times the 
size of the starburst nucleus, able to gen- 
erate X-rays. In this framework it is hard to 
understand the extended (more than 3 kpc 
long) X-ray emission of M82, a source that 
presents a more than 10 kpc long open chan- 
nel into the intergalactic medium, and that 
clearly is in the supergalactic wind stage. 

• Models with a single source predict also a 
strong Hq, omission arising from the lateral 
walls of the disrupted suporshell, continu- 
ously impacted by the UV radiation aris- 
ing from the central source (see for example 
Suchkov et al. 1994, Tenorio-Tagle & Muiioz 
Tuhon 1997, 1998). The emission in such a 
case should be strongly limb brightened. 

• Most models (perhaps with the only excep- 
tion of Tenorio-Tagle & Muiioz Tuhon 1997, 
1998, which accounts for the infall of mat- 
ter into the central starburst) also end up 
with a remnant that presents a wide open 
waist along the galaxy plane (see figures in 
Tomisaka & Ikeuchi, 1988; Suchkov etal.. 



1994). This could measure several kpc in 
radius and is very different to the 150 pc ra- 
dius estimated for the central perturbed area 
of M82. 

The Ha emission of M82 is not limb brightened 

and arises from distinct filaments that emanate 
from the base of the outflow. This last point is 
also relevant as many models attempt to explain 
the filamentary structure with shell instabilities 
or with matter entraining the supershell at high 
distances from the galaxy plane (see Shuchkov et 
al. 1994). 

In the following sections we derive the prop- 
erties of supergalactic winds powered by multiple 
superstar clusters and stress the main differences 
in the resultant structure with respect to models 
that assume a single source. 

2.1. The high metallicity outflow of super- 
galactic winds 

Weather one considers one or multiple sources, 
the amount of metals expected from massive star- 
bursts and thus to be found within SGWs, de- 
pends strongly on the assumed stellar evolution 
models. Two extreme possibilities were investi- 
gated by Silich et al. (2001) taking into consid- 
eration models with and without stellar winds by 
Maedcr (1992), Woosley et al. (1993) and Thiele- 
mann et al. (1993) (see also Pilyugin & Edmunds 
1996). The results indicate, first of all, that al- 
most 40% of the mass gone into stars is violently 
returned to the surrounding medium by means 
of winds and supernovae. Of this, about 4% is 
in oxygen in the case of models without winds, 
and 1% for models that include winds. In either 
case, if one assumes that the metals mix efficiently 
with the stellar hydrogen envelopes of the progen- 
itors, within the radius that encompasses the stel- 
lar cluster, then the metallicity of the superwind 
can be calculated. As shown by Silich et al. (2001) 
for the case of superbubbles, this reaches superso- 
lar values in all considered cases, showing a max- 
imum within the first 7-8 Myr of the evolution. 
Such values have recently been confirmed through 
X-ray observations by Martin et al. (2002) for 
NGC 1569. 

The large metallicities expected from massive 
star clusters are to have a strong impact on the 
cooling properties of the outflow and thus also 
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on the observational properties of freely expand- 
ing superwinds. Figure 1 shows the run of the 
expected metallicity of the matter ejected by a 
massive burst of stellar formation, using oxygen 
as tracer and stellar evolution models with winds, 
as a function of time. The estimate is for a co- 
eval starburst powered during 40 Myr of evolution 
(until the last 8 Mq star explodes as supernova 
in a coeval starburst model). The outflow, before 
the supernova era. will display a metallicity simi- 
lar to that of the gas cloud out of which the star- 
burst formed. Thus the gas ejected first presents 
the metallicity here assumed for the host galaxy 
{ZjsM = 0.1 Zq). The supernova products how- 
ever, rapidly enhance the metallicity of the out- 
flow, reaching values well above Zq for at least 
20 Myr of the evolution. Afterwards, once the 
yield is reached, the outflow steadily approaches 
the metallicity values of the host galaxy. Starburst 
models, with single or multiple energy sources, are 
to account for the high metallicities emanating 
from the massive centers of star formation, to de- 
rive in a consistent manner the impact of radiative 
cooling in the resultant outflows. 

2.2. Boundary and initial conditions 

Several two dimensional calculations using as 
initial condition CC85 adiabatic flows have been 
performed with the explicit Eulerian finite differ- 
ence code described by Tenorio-Tagle & Muiioz- 
Tunon (1997, 1998). This has been adapted to al- 
low for the continuous injection of multiple winds 
(see below). 

Here we consider the winds from several iden- 
tical SSCs, each with a mechanical energy depo- 
sition rate equal to 10**^ erg s~^. The energy is 
dumped at every time step within the central 5 
pc of each of the sources following the adiabatic 
solution of Chevalier & Clegg (1985). The separa- 
tion between the sources is clearly arbitrary. We 
have considered two different configurations. The 
first one with three SSCs placed at 60 pc and 90 pc 
from the central one, all of them sitting on a plane. 
A second configuration considers only two sources 
sitting at 30 pc and 60 pc away from the symme- 
try axis. Several calculations were made to reasure 
that the spatial resolution used, led to a conver- 
gent solution. All calculations here presented were 
made with the same numerical resolution of half a 
pc and all of them with an open boundary along 



the grid outer edge. The time dependent calcula- 
tions do not consider thermal conductivity but do 
account for radiative cooling, with a cooling law 
(Raymond et al. 1976) scaled to the metallicity as- 
sumed for every case. Here we present the results 
of different cases for which the assumed metallic- 
ity of the winds was set equal to 3Zq and IOZ0, 
justified by the high metallicity outflows expected 
from massive bursts of star formation (see Figure 
l). 

In all cases it is assumed that at the heart of 
each SSC, within the region that encompasses each 
of the rc;c;c!ntly formed stellar clusters (Rsb), the 
matter ejected by strong stellar winds and super- 
nova explosions is fully thermalized (CC85, see 
also Canto et al. 2000 and Raga et al. 2001). 
This generates the large overpressure responsible 
for the mechanical luminosity associated to each 
of the super clusters. Within each star cluster re- 
gion, the mean total mechanical energy Lsb and 
mass MsB deposition rates, control, together with 
the actual size of the star forming region Rsb, 
the properties of the resultant outflow. The total 
mass and energy deposition rates define the cen- 
tral temperature Tsb and thus the sound speed 
csB at the cluster boundary. As shown by CC85 
at the boundary r = Rsb, the flow starts expand- 
ing with its own sound speed. There is however 
a rapid evolution and as matter streams away it 
is immediately accelerated by the steep pressure 
gradients and rapidly reaches its terminal velocity 
(Vcx) ~ 2csb). This is due to a fast conversion 
of thermal energy into kinetic energy of the re- 
sultant winds. In this way, as the winds expand, 
their density, temperature and thermal pressure 
will drop as r~^, r~^/^ and r~^^^^, respectively 
(sec CC85). The flow is then exposed to suffer 
multiple interactions with neighboring winds and 
as shown in Paper I, it is also exposed to radiative 
cooling. For the former, the issue is the separation 
between neighboring sources and for the latter the 
local values of density, temperature and metallic- 
ity. Radiative cooling would preferably impact the 
more powerful and more compact sources, leading 
to cold (T ~ 10^ K) highly supersonic streams (see 
Paper I). 

2.3. The structure of supergalactic winds 

Figure 2 compares the initial stages of cases 1 
and 2, that consider three equally powerful {Lsb 
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= 10^^ erg s~^) superstar clusters sitting at 0, 60 
and 90 pc from the symmetry axis. All of them 
with an Rsb = 5 pc, produce almost immediately 
a stream with a terminal velocity equal to 1000 km 
s~^. The only difference between the two cases is 
the assumed metallicity set equal to 3Zq in case 1 
(upper panels) and lOZ© in case 2 (lower panels). 
At t = yr the three clusters are embedded in a 
uniform low density (p = 10~^^ g cm~^) medium. 
Thus our calculations do not address the devel- 
opment of a superbubble, nor the phenomenon 
of breakout from a galaxy disk or the halo, into 
the IGM. The initial condition assumes that prior 
events have evacuated the region surrounding the 
superstar clusters, and we center our attention on 
the interaction of the supersonic outflows. 

Figure 2 shows the resultant initial stages of 
cases 1 and 2. The various panels display the run 
of density and velocity (left panels) and that of 
temperature (next four panels) in four tempera- 
ture ranges: The regime of H recombination 10* 
K - 10^ K, followed by two regimes of soft X-ray 
emission 10^ K - 10^ K, and 10'^ K - 10^ K and 
the hard X-ray emitting gas with temperatures be- 
tween lO'' K - 10^ K. 

The crowding of the isocontours in the figures 
indicate steep gradient both in density or in tem- 
perature and velocity, and thus indicate the pres- 
ence of shocks and of rapid cooling zones. In the 
temperature plots one can determine the distance 
(~ 30 pc in case 1 and 15 pc in case 2) within 
the diverging outflows emanating from each of the 
superstar clusters, at which strong radiative cool- 
ing (in agreement with our analytical estimates in 
Paper I) becomes important in the two cases. 

The interaction of neighboring supersonic 
winds causes the immediate formation of their 
respective reverse shocks, and of a high pressure 
region right behind them. The pressure (and tem- 
perature) reaches its largest values at the base of 
the interaction plane, exactly where the reverse 
shocks are perpendicular to the incoming streams. 
The high pressure gas then streams into lower 
pressure regions, defining together with radiative 
cooling, how broad or narrow the high pressure 
zones, behind the reverse shocks, are going to be. 
Radiative cooling occurs in every parcel of gas at 
the rate prescribed by its local density, tempera- 
ture and metallicity. If cooling is avoided at least 
partially or temporarily, as in the first case, the 



high pressure region between the reverse shocks 
would drive them against the incoming streams, 
and very shortly they would acquire a oblique 
standing stable configuration to be retained for as 
long as the winds continue to interact (see Figure 
2, upper panels). 

This also happens if cooling is fast enough, the 
oblique reverse shocks rapidly acquire a standing 
location, however in these cases, the loss of tem- 
perature behind the shocks is compensated by gas 
condensation, leading, as in the second case (see 
Figure 2, bottom panels), to narrow, dense and 
cold filaments. The drastic drop in temperature 
occurs near the base of the outflow, where the gas 
density is large and radiative cooling is exacer- 
bated. The dense structures are then launched 
at considerable speeds (~ several hundreds of km 
s~^) from zones near the plane of the galaxy (see 
Figure 2 lower panels). These dense and cold 
structures are easy target to the UV radiation pro- 
duced by the superstar clusters and thus upon 
cooling and recombination are likely to become 
photoionized. Note however that as the free winds 
continue to strike upon these structures, even at 
large distances from their origin, the resultant cold 
filaments give the appearance of being enveloped 
by soft X-ray emitting streams. 

All of these shocks are largely oblique to the 
incoming streams and thus lead to two major ef- 
fects: a) partial thermalization and b) coUimation 
of the outfiow. These effects result from the fact 
that only the component of the original isotropic 
outflow velocity perpendicular to the shocks is 
thermalized, while the parallel component is fully 
transmitted and thus causes the deflection of the 
outflow towards the shocks. This leads both, to 
an efficient coUimation of the outflow in a general 
direction perpendicular to the plane of the galaxy, 
and to a substantial soft X-ray emission associated 
with the dense filamentary structure, extending 
up to large distances (kpc) from the plane of the 
galaxy. In the figures one can clearly appreciate 
that the oblique shocks, confronting the originally 
diverging flows, lead to distinct regions where the 
gas acquires very different temperatures, and thus 
that will radiate in different energy bands. 

Figures 3 and 4 show the time sequence of cases 
1 (with Z = 3 Zq) and 2 (with Z = 10 Z©), respec- 
tively, until they reach dimensions of one kpc, to- 
gether with the final temperature structure split- 
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ted into the four temperature regimes considered 
earlier (last four panels). 

The stream of gas behind the reverse shocks 
leads eventually to the establishment of crossing 
shocks at the tips of the oblique structures (see 
Figure 3), which are to become also stationary as 
the flow is effectively channeled into the IGM. 

As in the case of colliding stellar winds from 
binary systems (see Stevens, Blondin & Pollock 
1992) a variety of dynamical instabilities are found 
to dominate the shocked region, particularly when 
strong radiative cooling sets in (see Figure 4). 
These lead to the various loops and twists along 
the dense filamentary structures, which neverthe- 
less do not impede that the outflow reaches large 
distances away from the galaxy plane. The loops 
and twists along these structures promote also a 
larger cross-section to the incoming free wind and 
partly thermalized wind and thus lead to regions 
of enhanced soft X-ray emission clearly associated 
to the twisted Ha filaments (see Figure 4). 

Figure 5 displays the results from a final case 
in which two superstar clusters sitting at 30 and 
60 pc away from the symmetry axis interact to 
shape the inner structure of a superwind. The 
calculation also assumes aZ = lOZ0. As in case 
2, elongated filaments result from the interaction 
of the high metallicity winds, channeling into the 
IGM most of the energy deposited by the SSCs. 
Note that, as in the preceding cases, about 50% 
or less of the energy deposited by the most outer 
SSC is lost in the radial direction, while the rest, 
as well as that deposited by other energy sources, 
is fully driven into the IGM. 

2.4. Self-collimated supergalactic winds 

The high degree of coUimation attained in our 
calculations, that composes a SGW from a col- 
lection of energetic neighboring superstar clusters, 
results from the simple fact of having placed the 
individual wind sources of equal strength, all of 
them, in a preferential plane. In this way, it be- 
comes irrelevant if they all sit on a flattened disk, 
or a ring. As long as they all sit near a pref- 
erential plane, the interaction of neighboring su- 
personic diverging flows will promote the multi- 
ple standing reverse oblique shocks and crossing 
shocks that will unavoidably lead to a remarkably 
eflicient self-coUimation. CoUimation that does 



not required of a torus or a thick disk of ISM. 
If all superstar clusters sit on a plane, only a frac- 
tion of the energy provided by the ones sitting at 
the most outer extremes of the cluster distribution 
will interact with the general ISM. However, most 
of the energy produced by the collection of sources 
will be rechanneled by the standing oblique shocks 
resulting from neighboring interactions, to com- 
pose a broad base supersonic jet, capable of self- 
collimation. The base of the outflow will then have 
dimensions similar to the flattened cluster distri- 
bution and as shown above, depending on the in- 
dividual energetics, proximity and metallicity, the 
general outflow is to generate a dense and cold fll- 
amentary structure as well as a kpc extended soft 
X-ray emitting region. A rich structure that could 
not arise if one assumes a free expanding wind that 
emanates from a single SSC. 

A wider jet structure may be generated in cases 
in which the population of SSCs do not have the 
same mass and thus equal mechanical energy input 
rates. Under such conditions the oblique shocks 
will present standing configurations more inclined 
over the less energetic clusters and this will lead 
to the fanning and broadening of the outflow, and 
to the inclination of the fllamentary pattern. 

3. Discussion 

From the starburst synthesis models (see Lei- 
therer & Heckman 1995) one knows that a super- 
star cluster with a total mass in stars (say M« 

10^ Mq) produces an almost constant ionizing 
photon flux {Fuv ^ 10^^ photons s"'^) during the 
flrst few (3-4) Myr and then abruptly, it begins to 
decrease as t~^ as the most massive stars begin 
to evolve away from the main sequence to eventu- 
ally end up as supernovae. This implies that after 
10 Myr of evolution, the ionizing flux would be 
more than two orders of magnitude smaller than 
its original value and thus the UV radiation will 
be unable to ionized the original HII region vol- 
ume, limiting to 10 Myr the duration of the HII 
region phase. On the other hand, the mechani- 
cal energy deposition from such a cluster leads to 
an almost constant value (^ 10"*° erg s"'^) over a 
much longer time span, as it includes the corre- 
lated supernovae from stars down to 8 Mq with 
an evolutionary time of 40 - 50 Myr. And thus 
the supernova phase is 4 or 5 times longer than 
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the HII region phase. 

Under such circumstances, if one considers a 
starburst nucleus composed by several SSCs gener- 
ated at different times, then the time span during 
which the isotropic winds from these may interact, 
the coherence phase, is limited to 40 Myr. Within 
this time, newly born clusters will have the ca- 
pability of causing the ionization of the structure 
produce by interacting winds that emanate even 
from cluster with an age in excess of 10 Myr. Dur- 
ing the coherence phase, some of the SSCs will also 
be producing highly metallic outflows (see Figure 
1), the interaction of which will lead to a filamen- 
tary wind structure. 

A comparison of the last calculated time of 
cases 1 and 2, each with 3 SSCs dumping 10^^ 
erg (see Figures 3 and 4), when the redirected 
outflow has reached dimensions of almost 1 kpc, 
allows for an estimate of the flUing factor occu- 
pied by gas at different temperatures. The hot (T 
> 10^ K) gas occupies almost 70% in case 1 and 
30% in case 2, of the total area. The warm gas 
(T < 10^ K) fills most of the remaining volume al- 
though the dense and cold enhancements evident 
in Figure 4 occupy about 40% of the superwind 
cross-sectional view. There is also a small vol- 
ume around the SSCs that present temperatures 
that will allow the gas to radiate in the hard X-ray 
regime (see last two plots in Figures 3 and 4). This 
numbers are to be compared with the results from 
the outflow produced by an equally energetic (3 x 
IQii gj.g g-i^ single cluster. Following Paper I, we 
have calculated the temperature distribution and 
thus the radius at which radiative cooling (assum- 
ing Z = 3Zq) sets in for a clusters with a radius 
of 95 pc (the size of the cluster distribution used 
in cases 1 and 2). The temperature of such a free 
streaming outflow plummets to 10'' K at 480 pc 
(instead of ^ 10 kpc obtained if one assmncs an 
adiabatic flow). In such a case, if one considers 
a similar volume to that displayed in Figures 3 
and 4, then the hot gas (T > 10^ K) filling factor 
will be 37.5% and the rest of the volume will be 
filled with gas capable of being re-ionized by the 
stellar photon flux. When comparing the results 
from single and multiple sources, it is central to 
notice the spatial distribution of the various re- 
sultant gaseous phases. Cases with multiple SSCs 
lead to the spatial co-existence of the X-ray (T > 
10^) and the dense and cold (T < 10^ K) emit- 



ting filaments (sec panels 5 and 6 in Figures 3 and 
4). In the case of a single source of energy the X- 
ray emitting gas is not in direct contact with the 
cooler medium along the outflow, i. c. the struc- 
ture of the outflow is concentric, with the X-rays 
emanating only from the most central regions. 

The coUimation caused by the various oblique 
and crossing shocks in the multiple source cases, 
that makes the superwind avoid the diverging out- 
flow inherent to isotropic single source cases, chan- 
nels almost flve times more energy within the com- 
putational area above considered. In case 1 the 
thermal and kinetic energy of the hot phase domi- 
nate with 26% and 57%, respectively, over the 17% 
kinetic energy found in the gas with T < 10^ K. 
These numbers are to be compared with the re- 
sults from the single source case described above 
that present within a similar computational vol- 
ume, a 68% and 27% as thermal and kinetic en- 
ergy of the hot gas, while only a 5% of the total 
appears as kinetic energy of the gas with T < 10^ 
K. 

The origin of the X-rays in nuclear starburst re- 
gions and of the filamentary structure, as seen in 
Ha in M82, have been ascribed to features seen in 
models that consider various stages in the devel- 
opment of hot superbubbles powered by a central 
starburst. Models that present a single reverse 
shock, an outer super shell and thus models that 
have little to do with a free expanding SGW. Note 
that if the wind of M82 has reached the Ha cap 
at 10 kpc from the galaxy disk, and is expand- 
ing with say. 1000 km s^^, then the free stream- 
ing outflow started at least 10 Myr ago. During 
that time the base of the outflow has managed to 
preserve a comparatively small dimension (radius 
~ 150 pc) implying a very efflcient channeling of 
the deposited energy in a direction away from the 
disk of the galaxy. There is also a growing evi- 
dence of large-scale features in starburst galaxies, 
caused by an important stellar energy input rate 
and the richness of structure within the ISM. How- 
ever, in most of the cases the evidence for a freely 
expanding supergalactic wind is only marginal. In 
the case of NGC 253, its extended X-ray emit- 
ting bubble is much smaller than the dimensions 
of the dusty galaxy halo found by Melo et al. 2002, 
implying therefore that it is still evolving along 
the superbubble phase. Similar conclusions were 
drawn by Martin et al. (2002) for NGC 1569: 
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"The X-ray color variations in the halo are incon- 
sistent with a free-streaming wind and probably 
reveal the location of shocks created by the inter- 
action of the wind with a gaseous halo" . Several 
more dwarf starburst galaxies were considered by 
Legrand et al. (2001) where the energetics inferred 
from the central clusters were compared with the 
limit for mass ejection derived by Silich & Tenorio- 
Tagle (2001). In all considered cases the bulk of 
the galaxy sample lie below the limit required to 
reach the galaxies outskirts. 

Prom the results of section 2, it is clear that the 
inner structure of super galactic winds strongly de- 
pends on the energy and mass deposition history. 
In particular we have shown that the richness of 
structure is largely enhanced when the presence of 
superstar clusters, their powerful winds, and pos- 
sible interaction within a single starburst nucleus, 
are taken into consideration. These considerations 
open a new set of possibilities. Issues such as 
the intensity of star formation in every superstar 
cluster, which defines their mechanical luminosity, 
their age which also impacts on the metallicity of 
the ejected matter, as well as the number of su- 
perstar clusters, their compactness, and their po- 
sition within a starburst nucleus, are all relevant 
new parameters that allow for the co-existence of 
X-rays and optical emitting features, even at large 
distances from the source of energy. 

From our results it is clear that a plethora of 
structure, both in X-rays and in the optical line 
regime, may originate from the hydrodynamical 
interaction of multiple winds. The interaction 
leads to multiple standing oblique (reverse) shocks 
and crossing shocks able to coUimatc the outflow 
away from the plane of the galaxy. In our two di- 
mensional simulations, these are surfaces that be- 
come oblique to the incoming streams and thus 
evolve into oblique shocks that thermalize only 
partly the kinetic energy of the winds causing a 
substantial X-ray emission at large distances away 
from the galaxy plane. Surfaces that at the same 
time act as collimators, redirecting the winds in a 
direction perpendicular to the plane occupied by 
the collection of SSCs. Radiative cooling behind 
the oblique shocks leads, as soon as it sets in, to 
condensation of the shocked gas, and thus to the 
natural development of a network of filaments that 
forms near the base of the outflow, and stream 
away from the plane of the galaxy to reach kpc 



scales. Under many circumstances these filaments 
develop right at the base of the outflow and for all 
cases the prediction is that they are highly metal- 
lic. Note that the speed with which the calculated 
filaments raise above the galaxy plane is ~ 600 
km s~^, a value in excellent agreement with the 
measured deprojected velocities of the filaments 
in M82 (ShopbeU & Bland-Hawthorn, 1998). Hy- 
drodynamic instabilities play also a major role on 
the filamentary structure. Non-linear thin shell in- 
stabilities as studied by Vishniac (1994) as well as 
Kelvin Helmholtz instabilities broaden, twist and 
generally shape the filaments as these stream up- 
wards and reach kpc scales. The broadening of 
the filaments causes their interaction with the free 
winds, thermalizing further the rapid stream while 
leading to the development of soft X-ray emitting 
zones that envelope the densest structures. 

The surfaces that develop at the plane of inter- 
action between two wind sources are in our two 
dimensional approach depicted as vertical struc- 
tures. From these, the only real vertical structure 
is the filament that forms along the symmetry axis 
in our last case. This results from the convergence 
of multiple winds arising from superstar clusters 
sitting on a ring 30 pc away from the symmetry 
axis. It is indeed necessary to perform our calcula- 
tions in three dimensions to see the final outcome. 
Note however that the lateral side of all interaction 
planes will be launched into the highest pressure 
regions i.e. close to the base of the interaction, 
and thus it is very likely that they would be de- 
stroyed by the collision with other similar struc- 
tures arising from other interaction planes. The 
final outcome is thus expected to be very similar 
to that depicted by our two dimensional simula- 
tions. Nevertheless three dimensional simiilations 
are now underway and will also consider a variety 
of stellar masses and ages of the superstar clusters, 
as well as diflFerent locations and numbers within 
a starburst region. 
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Figure captions. 

Fig. 1. — The metallicity of the ejected matter. 
The metallicity of freely streaming outflows pro- 
duced by massive burst of star formation (in so- 
lar units) is plotted against the evolutionary time. 
The outflow emanates from a coeval starburst able 
to thermalize and mix all the newly processed met- 
als with the stellar hydrogen envelopes of the pro- 
genitors. The curve is derived for stellar evolution 
models with winds (see Silich et al. 2001) using 
oxygen as tracer. The metallicity of the outflow 
reaches super solar values during most of the evo- 
lution and is strongly reduced down to the original 
assumed ISM values, once the production of oxy- 
gen reaches its yield. 

Fig. 2. Two dimensional supcrwinds. The var- 
ious panels in every row represent cross-sectional 
cuts along the computational grid showing: iso- 
density contours with a separation A log p = 0.1 
and the velocity field for which the longest arrow 
represents 10^ km s~^. The following four panels 
display isotcmpcraturc contours, within the range 
10^ K - 10^ K, IQS K - 10^ K, 10^ K - 10^ K and 
lO'' K - 10^ K, respectively. Each of the super- 
winds has a power of 10'*"'^ erg and a radius of 
5 pc. The evolution of each wind starts from the 
adiabatic solution of CC85. The distance between 
consecutive tick mark = 25 pc in all figures. Thus, 
the size of the plots is 100 pc x 250 pc. The evolu- 
tion time for the two models is: 1.62 x 10^ yr and 
1.79 X 10^ yr. The assumed metallicitics were Z 
= 3 Z© (upper panels), and 10 Zq (lower panels). 

Fig. 3.— The same as Figure 2 for Z = 3Z0. The 

evolutionary time in the four first panels is 1.62 
X 10^ yr, 4.17 x 10^ yr, 9.4 x 10^ yr and 1.25 x 
10® yr. The last four panels show the temperature 
distribution, as in Figure 2, for the last calculated 
model. The size of the plots displays the whole 
computational grid: 100 pc x 1 kpc. 

Fig. 4. The same as Figure 3 for Z — WZq. The 
evolutionary times of the first four panels is 1.79 
X 10^ yr, 4.82 x lO"^ yr, 1.05 x 10*^ yr and 1.39 x 
10^ yr, respectively. 

Fig. 5. — The same as Figure 3 for Z = lOZ©. 
The calculation considers only two superstar clus- 
ters far from the symmetry axis. The evolutionary 



times of the first four panels is 1.8 x 10^ yr, 4.94 
X 10^ yr, 1.05 x 10® yr and 1.4 x 10® yr, respec- 
tively. 
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